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From the experimental data, collected at HERA ep-collider ||], one can conclude that the charm 
production model proposed in article |^ (below denoted as BKL), allows to describe inclusive photopro- 
duction and deep inelastic production of £'*^(2010)-mesons (below denoted as iD*-mesons), as well as 
inclusive photoproduction Z?s-mesons, with a good accuracy. 

Recently new data on Z?* -meson diffractive photoproduction have been presented by the ZEUS col- 
laboration [Q. In this connection we try to describe these new data in the framework of the BKL model. 

Let us remind the general features of this model. In the BKL approach one needs to produce pertur- 
batively c- and d-quarks which softly form Z?*-meson. These perturbatively produced c- and d-quark are 
valence ones for the meson. The soft hadronization process of the cJ-pair color singlet state is described 
by the average value of the operator: 

^^(1)^ = 12^7^ {'^"^ + %f^) {D*{j>D.)\{c^,d){d^.c)\D*(j>D^)), (1) 

where p^. is Z?*-meson momentum, and Md* is the meson mass. In the framework of the nonrelativistic 
potential model the average value of the operator correspond with the squared wave function in the origin: 
{0[i))\nr — l^(0)p. The hadronization of the color octet state is described by the average value of the 
analogous operator: 

= {^^'^ + %|^) {D*{pD>)\{-c^,\''d){d^,\'c)\D*{pD>)) ^. (2) 

It is worth to mention that the BKL model is based on the partonic concept of the hadronic structure. 
Indeed, in the framework of the partonic model for the system of the infinite momentum the valence quark 
structure functions are determined as follows: 

fcix,PT) = fc{x,PT) - fc{x,PT), 

ni^'PT) = fdi^^Pr) - fd{x,pT), 

where px is a transverse momentum of a parton inside the hadron, and x is the fraction of the hadronic 
momentum taken away by the parton. The average of momentum fractions taken away by the valence 
quarks are described by the following equations: 

{xD = / (fpxdx X ■ f^{x,pT) 



{x}) = / d^prdx X ■ fj{x,pT) « 

where (x^) + (xj) ~ 1, and A is the energy of the quark coupling into the meson. In the framework of 
the BKL model we neglect the quark velocity difference and assume quark velocities to be equal to each 
other: Vc = vj. The effective mass of the light quark in such approach plays the role of the infrared 
cut. We take this mass equal to A. Thus we obtain the following equations: 

«) = ^ 



Now it is clear that the BKL model is an extention of the parton model for the case of final hadrons 
in the framework of the valence quark approximation. 

This fact distinguishes the BKL approach from perturbative calculations ||, ^ based on the frag- 
mentation model of the hadronisation. In the framework of the fragmentation model of the hadronization 
it is supposed that the perturbatively produced single c-quark becomes a meson at large distances. This 
meson gets a fraction z of the c-quark transverse momentum kx with the probability determined by the 
fragmentation function Dc^d*{z^ pi): 



dzdpfj^* dkr 



kT = ^ 



2 



where Dc— >d* (-2, m) is normalized to the probabihty of c-quark to become a Z?*-meson w{c D*), 
measured in the e+e^-annihilation [|| {w{c D*) = 0.22 ± 0.014 ± 0.014). /i is the scale at which the 
perturbative partonic cross section of the cc-pair production dacc/dkr is calculated. 

It is clear that in the fragmentation approach one can not take into account the possibility for the c- 
quark to hadronize via the interaction with the quark sea of the initial hadron (recombination mechanism). 
That is why one needs to account this opportunity for c-quark to become a Z)*-meson in the frame work 
of some additional model. 

In the framework of the BKL approach both fragmentation and recombination mechanisms are ac- 
counted naturally in the calculations. It is worth to mention that both the fragmetation mechanism and 
the recombination one are desctibed by the same set of the diagrams. 

The fragmentation mechanism dominates at the large transverse momentum of the Z)*-meson in 
accordance with the factorization theorem. The main contribution at small transverse momenta is due 
to the recombination mechanism, i.e. due to the fusion of c-quark and a light quark from the sea of the 
initial hadron. It is worth to mention that the recombination contribution mechanism corresponds with 
higher twist cotribution to the transverse momentum distribution. 

As it was shown in many articles, devoted to the fragmentation model, the particular form of the frag- 
mentation function affects calculation results insignificantly. In the majority of the articles the Peterson 
B fragmentation function is used: 



where N is normalization factor, and e is a free phenomenological parameter dependent on the scale /i. 
The results would not change crucially if the parameterization of Kartvelishvili-Likhoded-Petrov [0 is 
used: 

D{z) ^ Nz'^'-'il- zy-'^'', (8) 

where Uc = —3, ad — 1/2 and 7 = 3/2. 

In the BKL model at large tranverse momenta (p^ > 20 ) the D*-meson cross-section can be 
expressed by formula (H), if the following pertubatively motivated form for fragmentation function is 



used [0: 



2(1 - r)(4 - r + 2r'^)z^ + (1 - rfi^ - 2r + 2r'^)z% (9) 
where r =^ md/^nid + rric) and 

(0^^^) = (0(i)) + i(0(8)). (10) 

Meson production in e^e^-annihilation is caused by fragmentation mechanism only. This fact allows 
to define {O^-^^) using w{c ^ D*) and quark masses: 

w{c^D)^ / Dc^D'{z)dz = q •^(^), (11) 

m3 







where J(r) is determined in ||]. The value of w{c D*) is known from the experiment. 
So one can determine the value of {O*^^^) for fixed values of m^, nic and ^. Assuming 



nid 0.3 GeV, 
rric — 1.5 GeV and 



w{c^D*) = 0.22. 

one obtains 

(O^f^imD*)) = 0.25 GeV^ 
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The best description of the experimental data on the charm photoproduction and the deep inelastic 
charm production at HERA can be achieved if (0(8))/(0(i)} = 1.3. 

The mentioned values of the parameters were taken both for the presented calculations and for the 
calculations of the iI'*-meson nondiffractive production cross section Q. 

BKL model can be used for all values of £'*-meson transverse momentum in the contrast to the 
fragmentation model which is valid only for large transverse momentum (our estimation shows that 
fragmentation model can be used for > 20 GeV). 

We use the parameter values (|l^) to calculate the cross-section of the I?*-meson diffractive production 
in the frame work of the BKL model. 

We choose one of the known form of the Pomeron flux parametrization : 



f ( i\ 111 



6.38e~^l*l +0.424e~3|* 



(13) 



where t is the squared transfer momentum in the proton vertex, and a;jp is the momentum fraction of 
the proton carried away by the Pomeron. We neglect the amplitude dependence on t in the calculations. 

We suppose that the Pomeron consists of the gluons only and used two types of the gluonic distribution 
inside the Pomeron G(/3): 

Rn(R\ - / ^^(1 " ^) ~ Pomeron; 

^^^^^ "16(1- - "soft" Pomeron, ^^^^ 

where (3 is the fraction of the Pomeron momentum carried off by a gluon. 

In Fig. 1 the calculation of the differential distributions of the Z?*-meson production has been performed 
for the kinematic region investigated by ZEUS Collaboration|@: 130 < < 280 GeV, < 1 GeV^, 
> 2 GeV, [77^ I < 1.5, 0.001 < xp < 0.018, where W is the invariant mass of the photon-proton 
system, is the photon virtuality and 77^ is the pseudorapidity of the D*-meson. The value of rj^ 
is determined by the angle 9 between the initial proton direction and the _D*-meson in the laboratory 
system as follows: rj^' = — ln(tg|). 

The cross sections calculated for this kinematic region in the frame work of BKL model have the 
following numerical values: 



^BKL 



0.77 ± 0.02 - "hard" Pomeron; 
0.56 ± 0.03 - "soft" Pomeron. 



From the transverse momentum distributions one can conclude, that the calculation with the soft 
Pomeron predicts more rapid decreasing of the cross section with increasing transverse momentum than 
the calculation with the hard Pomeron. 

From the distribution in D* pseudorapidity one can see that the model of the hard Pomeron predicts 
the cross section values of the D* meson production in the forward direction essentially larger than the 
ones predicted by the soft Pomeron model. The predictions of these models for the £)*-meson backward 
production in the considered kinematic region are practically the same. 

The hard Pomeron model predicts the distribution maximum at small AI^ values. The soft Pomeron 
model calculation yields the distribution maximum at Mx over 20 GeV. 

In the considered range oi xp the differential cross sections calculated with different models of 
Pomeron, show quite different behaviour: hard Pomeron decreases with xp while soft Pomeron increases 
instead. 

The ZEUS GoUaboration plans to continue the analysis of the diffractive photoproduction of the D* 
meson. That is why we have calculated the total and the differential cross sections for 0.001 < a;jp < 0.2 
(Fig. 2): 

BKL _/ 1-64 ±0.02 - "hard" Pomeron; 

5.45 ± 0.02 - "soft" Pomeron. 



a 



and for the 0.001 < xp < 0.1 (Fig. 3): 



^BKL 



1.51 ± 0.03 - "hard" Pomeron; 
3.56 ± 0.03 - "soft" Pomeron. 
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It is evident that for these ranges of xjp the difference between the soft Pomeron model and the hard 
Pomeron model is more essential than for the 0.001 < xp < 0.018. The soft Pomeron model leads to the 

larger value of the cross section than the hard Pomeron model. 

These researches have been partially supported by RFBR Grants 00-15-96645 and 99-02-16558. 
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Fig. 1. The BKL model predictions for the differential cross sections of the D* diffractive photo- 
production at HERA: a) p|? ; b) 77^ ; c) Mx] d) xp are calculated in the kinematic region: 
130<W < 280 GeV, < 1 GeV^, p!^' > 2 GeV, lr/^*| < 1.5, 0.001 < xp < 0.018. The 
solid histograms stand for the calculations with the hard Pomeron model, and the dashed 
ones stand for the calculations with the soft Pomeron one. 
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Fig. 2. The BKL model predictions for the differential cross section of the D* diffractive photopro- 
duction at HERA: a) ; h) r]^ ; c) Mx; d) xp are calculated in the kinematic region: 
130 <W < 280 GeV, < i CeV^, p^* > 2 GeV, \ri^'\ < 1.5, 0.001 < a;ip < 0.2. The 
solid curves stand for the calculations with the hard Pomeron model, and the dashed ones 
stand for the calculations with the soft Pomeron model. 
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Fig. 3. The BKL model predictions for the differential cross section of the D* diffractive pho- 
toproduction at HERA: a) ; b) 77^ ; c) Mx are calculated in the kinematic region: 
130 <W < 280 GeV, < 1 QeV^, p^* > 2 GeV, \7j^'\ < 1.5, 0.001 < a;ip < 0.1. The 
solid curves stand for the calculations with the hard Pomeron model, and the dashed ones 
stand for the calculations with the soft Pomeron model. 
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